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Abstract 

A search for the decay B± — > a^7r° with the Oq~ decaying to an r/ and a 7r + was carried out at the 
Stanford Linear Accelerator Center using the BABAR detector coupled with the PEP-II collider. 
The analysis used a data sample comprised of approximately 252 million BB pairs collected at the 
T(4S) resonance. No signal was observed and a 90% confidence level upper limit on the branching 
fraction was set at 1.32 x 10 -6 . 
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1 INTRODUCTION 



The quark content of the ao mesons is a subject of debate. It has been conjectured that they may 
not be simple qq states but may have a more complex nature such as having a KK component, 
being a glueball, or being a mixture of two and four quark states. According to Delepine et al. [2], 
a measurement of the branching fraction (BF) of B + — > clqTt may help in resolving this problem, 
as the predicted branching fraction for the four quark model is expected to be up to an order of 
magnitude less than the two quark model. The predicted branching fraction for this is already 
small, of order 2 x 1CP 7 . 

Assuming the two quark model, the proposed Feynman diagrams for the dominant tree level 
decays are given in Fig. ^ The low expected branching fraction is explained by the fact that the 
colour suppressed diagram (b) is expected to dominate over the colour allowed diagram (a) since 
the colour allowed case is also doubly suppressed by G-parity and vector current conservation 3 . 
Due to this it is also not possible to use isospin arguments to relate the branching fraction for this 
mode to others in the same final state Dalitz plane, such as B + — > aQ-7r + . 




(a) (b) 



Figure 1: The proposed Feynman diagrams for the process B + — ► ci^ii® with (a) external (color allowed) 
production of the aj and (b) internal (color suppressed) production of the &^ . 

The aim of the analysis described in this paper is to better constrain the existing models by 
measuring the BF for the decay B + — > clqTT . The analysis focuses solely on the case where the o^j" 
decays to an r\ and a ir + , with the r\ decaying to two photons. This sub-mode accounts for ~ 40% 
of all r/ decays. An unbinned extended maximum likelihood fit (ML) method is used to extract 
the B + — > clqIt event yield. Throughout this document the conjugate decay mode B~ — > tt° is 
implied and ao refers to the ao(980) unless stated otherwise. 

2 THE BABAR DETECTOR & DATASET 

The analysis was performed using data collected with the BABAR detector |I] at the PEP-II 
asymmetric-energy e + e~ collider 5^ at the Stanford Linear Accelerator Center. Charged parti- 
cles are detected and their momenta measured with a 5-layer double sided silicon vertex tracker 
(SVT) and a 40-layer drift chamber (DCH) inside a 1.5 T solenoidal magnet. A quartz bar ring- 
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imaging Cerenkov detector (DIRC) complements the dE/dx measurement in the DCH for the 
identification of charged particles. Energies of neutral particles are measured by an electromag- 
netic calorimeter (EMC) composed of 6, 580 CsI(Tl) crystals, and the instrumented magnetic flux 
return (IFR) is used to identify muons. 

The data sample consists of 229.8 fb -1 collected at the T(4S) resonance. The number of BB 
pairs used in the analysis totals (252.2 ± 2.8) x 10 6 |0j. 

3 RECONSTRUCTION & SELECTION 

The analysis was based on the reconstruction of final states consisting of four photons and one tt + . 
All charged particle tracks were taken as tt + candidates unless they were identified as kaons by a 
particle identification algorithm based on DCH and DIRC information. Tracks were required to 
have momentum < lOGeV/c with transverse momentum > O.lGeV/c. They were also required to 
be composed from > 12 DCH hits and with a distance of closest approach to the interaction region 
within 1.5 cm in the transverse direction and 10 cm along the axis of the detector. 

In order to construct rj candidates, all pairs of photons in each event were combined. The 
individual photons were required to have energy between 0.05 and 10 GeV. The invariant mass of 
each candidate pair was then required to be between 0.515 and 0.569 GeV/ c 2 . This required mass 
range was optimised by maximising the ratio of signal to the square root of the background, S/V~B, 
using signal and background Geant4 Monte Carlo (MC) samples. 

The r] candidates were refitted to constrain their mass to the known value [7j and were then 
combined with the 7r + candidates to form candidates, which were required to have a mass 
between 0.8 and 1.2 GeV/ c 2 . The ciq mass selection was left relatively loose as this variable was 
used in the ML fit. 

Pairs of photons in the event were also used to form 7r° candidates. The photons were required 
to have energy between 0.03 and 10 GeV. The tt° candidates were required to be consistent with 
the mass of a ir°, i.e. in the range 0.115 to 0.150 GeV/ c 2 . The ir° candidates were also refitted to 
constrain their mass to the known value [Jj. 

B meson candidates were formed by combining all a,Q and 7r° candidates in the event. The B can- 
didates were described kinematically using the energy substituted mass tties = [(^s + Po-Pb) 2 / Eq — 
\pb\ 2 ]z and the energy difference AE = E* B — ^\/s. Here, ps and Eb are the momentum and 
energy of the B candidate, po and Eq are the momentum and energy of the initial (e + ,e~) state 
and the * indicates this quantity is expressed in the centre-of-mass system, ^/s is the CM energy. 
Both m-Es and AE were included in the ML fit with the requirements, 5.20 < ities < 5.29 GeV/c 2 
and \ AE\ < 0.35 GeV. Since event reconstruction is an imperfect process more then one candidate 
was produced for each B. The average number of B candidates per signal event was observed in 
signal MC to be 1.38. Further selection based on the quality of these candidates was not performed 
and all were included in the ML fit sample. The fit was designed, as explained and verified in 
section 5, to only pick up one candidate per event. 

4 BACKGROUNDS 

The major background in this analysis comes from random combinations of particles in continuum 
e + e~ — > qq events (q = u,d,s,c). These were primarily rejected by placing a requirement on the 
cosine of the angle between the thrust axis of the reconstructed B candidate and the thrust axis of 
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Table 1: The charmless modes identified as potentially contributing significantly to the background. The 
nominal contributions used in the fit are given. 



Decay Mode 


Expected Yield (B candidates) 


Charm B decays 




Charged b — > c 


322 ± 32 


Neutral b — ► c 


200 ± 20 


Charmless B decays 




B± -> 


78 ± 10 


B ± -> a±^° 

1 


58 ± 13 


B± -> p±r; 

B± -> ^(1450)7? 

B° -» 7r°7r° 


37±? 
25 ±25 
19 ±4 


Inclusive B — > X s 7 


18±S 


B° 77^° 


14 ±14 


B° - a V 


12 ±12 


B± -» flOpi 


6±6 


B± -> aJ(1450)7T° 


5±5 


B — > 71^ 7771"° (non-resonant) 


2±2 


B ± — > 7r ± 7r°7r° (non-resonant) 


2±2 



the rest of the event, | cos(0tb)\] the thrust was calculated in the CM frame. The distribution of 
this variable peaks sharply near 1.0 for the two jet-like combinations from continuum qq pairs and 
is approximately uniform for B meson pairs; | cos(^tb)| was required to satisfy | cos(^t_b)| < 0.594. 
This requirement was chosen by again maximising S/yB using signal and background MC samples. 

In order to further separate signal from continuum a Fisher discriminant J- jS] was used in the 
ML fit. The discriminant was a weighted linear combination of the absolute value of the cosine of 
the angle between the direction of the B candidate momentum and the beam axis, | cos(^t_b)| as 
defined above, and the Lq and L2 Legendre polynomial projections of the energy flow of the event 
with respect to the B candidate thrust axis. The Lq and L2 quantities were formed by summing 
over all the neutral and charged particles in the event which were not used to form the B candidate. 
The Fisher discriminant was required to satisfy — 3 < T < 1. 

The BB background was split into charmed and charmless decay components, both of which 
were accounted for in the fit. The charmed component was modelled by removing the charmless 
component from the BB MC sample. The charmless contributions were selected by studying B 
decays in MC with identical or similar final states to the signal decay. Further contributions were 
identified from general charmless B decay MC samples which pass the selection criteria. These 
modes were modelled using exclusive MC samples. From the final state Dalitz plot we expect only 
1 significant peaking background (B^ — ► p ± (1450)r/) to overlap our signal. The other modes which 
contribute are mainly present through some degree of mis-reconstruction. All B background yields 
were held fixed at their expected values in the final fit. The expected contributions for these modes 
are listed in Tabic ^ 
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5 MAXIMUM LIKELIHOOD FIT 



The input variables to the extended ML fit were jtieSj AE, J 7 , and the aj" resonance mass. The 
probability density functions (PDF) used to model the data are discussed below, and were de- 
termined from MC simulation. The exception to this is the continuum component where most 
parameters were left free in the fit, thus enabling us to extract the parameters directly from the 
data. The only other floated quantities in the fit were the signal and continuum yields. 

5.1 Signal Model 

The peaking components in AE and ttjes were modelled with a Novosibirsk function [S] and two 
independent Gaussians respectively. The T shape was modelled using an asymmetric Gaussian 
and the ajj~ resonance mass with a Breit-Wigner. The resonance was modelled in signal MC with 
a mass of 0.98 GeV/c 2 and a width of 80 MeV/c 2 P3J|. 

The signal shape can be distorted by true signal candidates which have been mis-reconstructed 
and is referred to as self-crossfeed (SXF). The B candidates affected have more background- like 
distributions in the fit variables. If improperly accounted for these can result in a reduction of the 
discriminating power of the fit. 

In some signal events one true B candidate was seen as well as a number of SXF B candidates. 
As stated in Section 3 this gives an average of 1.38 B candidates per event. In order to remove 
any distortion which may result from SXF, therefore giving the purest possible signal shape, a 
model was created which combines the expected true signal shape (described above) with a model 
to describe the SXF. High purity samples of signal and SXF were used for this purpose. The 
true signal sample was obtained from MC by requiring that all generated daughter particles of the 
signal decay were correctly reconstructed. This method was not 100% efficient, leading to cross- 
contamination between the signal and SXF samples. In order to minimise this effect the signal and 
SXF samples were fitted iteratively to determine the PDF parameters. The SXF model was not 
explicitly included in the final fit since its purpose was to facilitate the removal of SXF distortion 
from the signal shape. The resulting PDFs, projected onto signal MC, are presented in figure [2j 

It was verified that there was negligible contamination of the signal component with SXF 
events and that these were instead absorbed by the continuum PDF. It was also verified that this 
method yields no greater than one good signal candidate per event, despite the inclusion of multiple 
candidates. The tests are outlined in Section 1531 

5.2 Background Model 

Any peaking contributions in toes an d AE in the background samples were modelled using Gaussian 
shapes. Where a slowly varying, combinatorial component exists it was modelled with low-order 
polynomials in AE and an ARGUS threshold function for toes- I n most charmless modes cor- 
relations were expected to exist between these variables. In order to model these correlations the 
above prescription was, where necessary, replaced with a two-dimensional non-parametric KEYS 
PDF 12 . T was modelled either with two Gaussians or an asymmetric Gaussian. Peaking com- 
ponents in the aj" mass variable were modelled with Gaussians or Breit-Wigners depending on 
mis-reconstruction. Any more slowly varying components were modelled with low-order polynomi- 
als. 
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AE (GeV) m ES (GeV/c 2 ) 




fisher discriminant a mass (GeV/c 2 ) 

Figure 2: Signal model PDFs projected onto signal MC. The total PDF is given in blue (solid line) with 
the true signal contribution in red (dashed) and the SXF in black (dash-dotted). Only the true signal shape 
is used in the final fit. 



5.3 Fit Validation 

In order to verify the consistency of the fitting procedure MC simulations were generated from 
the PDF model. By refitting to these datasets and assessing the shifts in the fitted signal yields 
any problems with the model could be detected. In the studies no significant deviation from 
the generated signal yield was found. In order to detect biases from any improperly modelled 
correlations among the variables the above study was repeated with embedded randomly selected 
events from the fully simulated signal and background MC samples. No significant biases from any 
of the background modes were found, nor from any correlations in signal MC. Any small biases 
which exist were used to estimate a systematic error contribution. Finally, in order to test the 
ability of the fit to separate true signal from SXF, SXF candidates were explicitly embedded from 
the signal MC and a check for shifts in the fitted signal yield was made. Once again the fitted 
signal yields were consistent with those expected. 

6 SYSTEMATIC STUDIES 

The systematic uncertainties fall into three categories: firstly the errors in the number of events 
extracted from the fit, secondly the errors in the efficiency, and lastly the errors in the BF result 
due to uncertainties in the total number of B mesons and in the daughter BFs. 
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6.1 Systematic errors derived from the fit 

Uncertainties resulting from the ML fit give additive systematic errors in the number of signal 
events derived from the fit. 

• Uncertainties in the PDF parameters used in the fit 

This uncertainty lies in the statistical errors in the fitted parameters of the PDFs which make 
up the model. These were due to the limited amount of available MC simulated events upon 
which to model any given PDF. Each parameter in a PDF was considered independent of the 
others, thus neglecting any correlations which may exist. The PDFs considered were signal, 
BB events with charmless modes removed, and continuum. The charmless contribution was 
expected to be small and was neglected. With the exception of the signal ao mass Breit- 
Wigner width, all of the parameters were varied within their errors and the final fit repeated. 
In all cases the resulting shifts in the fitted signal yield were taken as the systematic error 
contribution. For the signal ao mass Breit-Wigner width lower and upper bounds of 50 MeV/c 2 
and 100MeV/c 2 were assumed. This was to account for the uncertainty in the modelled ao 
lineshape. 

• Uncertainties in charmless B decay contributions 

The uncertainty here lies in the normalisation of the separate charmless backgrounds due to 
errors in the BFs used to estimate their contribution. These values are listed in Tabled The 
systematic error contribution was calculated as the shift in the fitted signal yield in the final 
fit to data where the contributions from each specific mode were varied within the errors of 
their BFs. For the cases where only an upper limit exists, the contributions were calculated 
assuming 50% of the limit as a central value. The systematic error was then estimated by 
varying this assumed central value by ±100%. 

• Uncertainties in the yields of the B — ► charm modes 

This uncertainty was ascertained by varying the yields of the B to charmed modes by ±10% 
as a conservative estimate. This yield uncertainty was dominated by the error in the efficiency 
for reconstructing these B modes and not the cross-section for their production. Once again 
the shifts in the signal yield were taken as systematic errors. 

• Uncertainties in the bias from the fit 

To estimate any potential fit bias, a simulation study was run embedding all B backgrounds 
from the MC and generating continuum MC according to the values of the continuum shape 
parameters resulting from the final fit. Zero signal events were embedded or generated. The 
systematic error was estimated as the sum in quadrature of 50% of the fitted bias and its 
statistical error. 

6.2 Systematic errors in the efficiency 

There are a number of sources of systematic uncertainty which affect the efficiency and are thus 
applied as a multiplicative correction to the final result. 

• Tracking and neutrals efficiencies 

This source of uncertainty arises in tracking where a global per track systematic error of 0.5% 
was assigned based on results of dedicated studies. There was also an uncertainty in the 
efficiency for reconstructing neutral particles. A systematic error of 3% each for every 7r° and 
rj in the final state was assigned. 
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• Data/MC agreement in the | cos 9tb\ variable 

Due to the imperfect agreement between data and MC samples, the selection criteria for 
\ cos 6tb\ require the assignment of a systematic error. Control sample studies have shown 
that tighter selections incur larger errors. As such a conservative 5% error due to this selection 
was assigned. 

• Statistical error 

Finally, the error due to limited MC statistics in the efficiency had to be accounted for. This 
was simply the binomial error in selecting a given number of events from a larger sample. 

6.3 Systematic errors contributing to the branching fraction 

These systematic contributions are also multiplicative errors. 

• Total number of BB events 

The total number of BB events in the data set was estimated to be (252.2 ± 2.8) x 10 6 . The 
error was taken clS EL systematic. 

• Uncertainties in the daughter decay BFs 

The errors in B(a~Q — > t]tt + ) and B(r\ — I 77) were taken from the Particle Data Group [7]. 
The error in B(a^ — ► t]tt + ) was calculated by taking the ratio of the partial widths of the two 
dominant decay modes; — ► t]tt + and — > KK. The measured ratio is 0.183 ± 0.024 
which, assuming all other decay modes are negligible, gives B(clq — > r]ir + ) = 0.85 ±0.02. The 
value taken for B(rj -> 77) was 0.3943 ± 0.0026. 

6.4 Summary of Systematics 

The results of all of the studies to estimate systematic error are presented in Table El 
Table 2: Estimated systematic errors in the final fit result. 



Source of Uncertainty 


77 -> 77 


Additive (Events) 




Fit Parameters 


+ 7.7 
-4.4 


Charmless Yields 


+2.3 
-1.5 


Charm Yields 


+0.2 
-0.0 


Fit Bias 


±1.6 


Total Additive (Events) 


+8.2 
-4.9 


Multiplicative (%) 




Neutral efficiency 


±6.0 


Tracking efficiency 


±0.5 


\cos(8tb)\ Selection 


±5.0 


MC Statistics 


±0.9 


Number of BB Events 


±1.1 


Daughter ao Decay BF 


±2.0 


Daughter 77 Decay BF 


±0.7 


Total Multiplicative (%) 


±8.2 
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7 RESULTS 



The results from the fit to the full data set are given in Table The yield is found to be consistent 
with zero. 

Table 3: The results of the fit to the full data set and other values required for calculating the branching 
fraction. All B background yields were held fixed to the values listed in Table 1. The upper limit is shown 
first with only the statistical error and then with the total error. 



Required quantity/result 




Candidates to fit 
Signal Yield (events) 
Continuum Yield (candidates) 
ML Fit bias (events) 


36098 
-18±11 
35324±190 
2.55 


Accepted cff. and BFs 




e(%) 

B(V - 77) (%) 
B(a+ - m +) (%) 


16.18 
39.43 
84.5 


Results 




Branching Fraction (xl0~ 6 ) 
Upper Limit 90% C.L. (xl0~ 6 ) 
Upper Limit 90% C.L. (xlO" 6 ) 


-l.5±^(stat)± u -Z(syst) 
< 1.06 (statistical error only) 
< 1.32 (total error) 




90 p- 



5.21 5.22 5.23 5.24 5.25 5.26 5.27 5.28 5.29 



BaBar j 

preliminary - 




m ES (GeV/c 2 ) 



0.1 0.2 0.3 

AE (GeV) 



Figure 3: Projection plots in toes and AE comparing MC generated from the overall PDF (blue solid 
curve) with the data (black points). The samples have been background reduced by requiring the ratio of 
likelihoods C sig /[C sig + SC bkg ] to be > 0.6. 



The branching fraction is calculated using the following formula 

b= y ~ Yb (1) 

where Y is the signal event yield from the fit, Yb is the fit bias, e is the efficiency for the B decaying 
via the studied mode, N BB - is the number of produced BB mesons and HiBi is the product of the 
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Figure 4: Log likelihood scan curve showing the effect of convoluting in systematic errors. The red dashed 
curve shows the likelihood for statistical errors only while the blue solid curve shows the likelihood including 
systematic errors. 



daughter branching fractions. The decay rates for T(4S) — > B + B~ and T(4S) — > B°B° are assumed 
to be equal. 

Projection plots for the final fit results are presented in Figure |31 These compare a sample of 
MC generated from the overall PDF with the data. Both samples have had background reduced 
by requiring that the likelihood ratio C s i g /[C s i g + EC^g] for any event be > 0.6. 

The likelihood C was convoluted with a Gaussian to include the systematic errors. A scan was 
performed to give the variation in the likelihood for different input yields. The resulting curve 
in — 21n(£/£o) as a function of the branching fraction is shown in Fig. QJ Note how this plot is 
scaled relative to the value for the minimum negative log likelihood (Co). From this distribution 
the upper limit of 1.32 x 10 -6 was computed by integrating the associated likelihood function to 
find the yield which bounds 90% of the area under the curve. In doing this a starting branching 
fraction value of zero was assumed. 

This upper limit is reported in Table El with and without the systematic error to indicate the 
degree to which the limit is statistics dominated. The largest source of systematic error comes from 
the oo lineshape. 

8 SUMMARY 

A search for the decay B + — > ajj~7r° was carried out using 252 million BB pairs. An unbinned 
extended maximum likelihood fit was used to obtain the result. The fitted signal yield was consistent 
with zero and thus a 90% C.L. upper limit of 1.32 x 10 -6 was set on the branching fraction for this 
decay. Our sensitivity to this mode was therefore insufficient to make it possible to differentiate 
between the two and four quark ao structure models as defined by Delepine et al. 0. 
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